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Abstract 



The duality in the Chalker-Coddington network model is examined. We 
are able to write down a duality relation for the edge state transmission coef- 
ficient, but only for a specific symmetric Hall geometry. Looking for broader 
implication of the duality, we calculate the transmission coefficient T in terms 
of the conductivity axx and axy in the diffusive limit. The edge state scatter- 
ing problem is reduced to solving the diffusion equation with two boundary 
conditions (dy - ^d^U = and [d^ + ^^^^^dyU = 0. We find that the 
resistances in the geometry considered are not necessarily measures of the 
resistivity and p^x — WT^ (-^ — ^ — T) holds only when pxy is quantized. 
We conclude that duality alone is not sufficient to explain the experimental 
findings of Shahar et al and that Landauer-Buttiker argument does not render 
the additional condition, contrary to previous expectation. 



I. INTRODUCTION 

We consider the transport properties of a two-dimensional (2D) disordered strip con- 
nected to two disorder- free leads, as shown in figure 1. In the transverse direction, the 
system is restricted to a finite width of W by confinement potentials. The entire system is 
subject to a perpendicular magnetic field. In the leads the edge states are the only current 
carrying states. Therefore the problem can also be viewed as that of edge states scattering 
through a disordered region. From the viewpoint of the Landauer-Buttiker formula 
the DC transport properties are determined by the scattering matrix at the Fermi energy. It 
has been shown that the total transmission or the total back scattering of edge states leads 
to the quantization of Hall conductance ; from the bulk point of view the quantization 
requires the localization of the bulk states. The equivalence between the edge state descrip- 
tion and the bulk point of view has generated interesting discussions However, there 
has been no quantitative analysis as to how the edge state transmission coefficient relates to 
the bulk conductivity or resistivity. It has been shown that the longitudinal resistance 
Rxx of the above system is 

R h 



Rxx rj. ^2 ' 
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where T and R are transmission and reflection coeflicient of a single edge state. The foUowing 
relation 

_ L Rh 

Pxx — 

has also been used in a number of articles , where p^^ is the longitudinal resistivity, W 
and L is the width and length of the sample respectively. The latter relation is based on 
the assumption that the resistance is a measure of the resistivity. This assumption is not 
justifled. 

Our work is also motivated by the recent experimental and subsequent theoretical work 
by Shahar et al @,|l3] concerning the duality between phases on either sides of the quantum 



Hall transitions. Using the Chalker-Coddington network model [|ri[, we are able to obtain 
a previously speculated duality relation [|10| for the transmission coeflicient and the longi- 
tudinal resistance in the geometry specifled above. To decide whether or not the duality 
relation explains the experiment, which was done in a diflerent geometry, one again needs 
to resolve the relation between resistance and resistivity. 

This paper addresses two issues. First we are concerned with the duality in Chalker- 
Coddington model and how this duality manifests itself in resistance measurement. Secondly 
we render a microscopic calculation for the edge state transmission coeflicient in the diflusive 
limit in terms of the bulk parameters, the longitudinal and Hall conductivity a^^ and Gxy 
We flnd a non-trivial relation between resistance and resistivity. Although the calculation 
helps to put a restriction to the implication of the duality, it resolves an independent issue 
of its own. 

Section H serves as a brief review of the quantum linear response theory. There we 
set up the starting point for the microscopic calculations. We emphasize that for flnite- 
sized systems with phase coherence, all measured quantities are conductances or resistances 
and are in principle sensitive to the way the measurement is set up. We make contrast 
between the bulk current density and the Landauer-Buttiker scattering point of view point. 
In following both approaches in later calculations, we demonstrate their equivalence. 

In section HI we express the Hall and longitudinal resistance in terms of the trans- 
mission coeflicient using the Landauer-Buttiker formula. Our derivation makes explicit the 
involvement of probes and leads in the measurement. In section IV, we proceed to discuss 
the duality in the discrete Chalker-Coddington model, which leads to an inverse relation 
for the longitudinal resistance, but only for samples with reflection symmetry. We also look 
for the implication of the duality in the continuous limit. There the duality is between the 
phase with its bare Hall conductivity at a^y and the phase at n — a^^, where n is a integer. 
(Note: from now on we use /h as the unit for conductivity and h/e^ for resistivity.) The 
duality relation in the continuous limit translates to a relation between the renormalized 
conductivity: crxy{(^xy) + ^xyi"^ — (^xy) = ^ (here (Jxyi^^xy) denotes the renormalized axy as 
a function its bare value cr^^). We show that this duality relation alone is not enough to 
explain the experimental claim of 0. To do so an additional constraint between axx and 
axy is required. 

In section V we calculate the transmission coeflicient and the resistances in the pertur- 
bative limit {axx > !)• We show that the ideal leads aflect the outcome of the measurement 
by imposing a boundary condition on the electro-chemical potential and the Hall resistance 
is thus artiflcially flxed to a quantized value. We flnd that only under the condition that the 
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Hall resistivity is quantized the resistance is proportional to resistivity. The edge state scat- 
tering problem in the diffusive limit is reduced to a special boundary problem. Its analytical 
solution by conformal mapping is discussed in section VI. We conclude with a discussion 
on the missing connections between the duality found in experiment and that in existing 
models for the quantum Hall effects. 



II. TWO FORMS OF THE LINEAR RESPONSE THEORY 

The quantum mechanical linear response theory for non-interacting electron gas can be 
put in two forms. In one form one writes the local current density as a functional of the 
external field: 



(r) = I drV(r,r')^.(r'), (2.1) 



where cr^iy(r, r^) is the bilocal conductivity and can be expressed in terms of the single particle 
Green's function G^{E) = 1/{E — H ±ir]) where H is the single particle Hamiltonian. (For 
detailed form of a{r^ r') see [0). The relation between the current and the external field is 
generally non-local and in the presence of magnetic field (7^^^(r,r^) contains not only Fermi 
surface contribution, but also contributions from all energies below the Fermi energy. 



Due to a set of current conservation constraints |T2[ , the total current in a lead (say, the 



ith lead) to linear order depends only on the voltages in the leads, the V^s: 

I^ = g^JVJ. (2.2) 

where the conductance coefficients, the gijs^ are surface integrals of cr(r,r^) in the ith and 
jth lead: 



9^J 



= - 1 1 dS,'a{r,r')'dSj. (2.3) 



Although the off-Fermi-surface terms in a(r,r') contributes to the local current response, 
they give zero net contribution upon surface integral [|12|. As a result the gijs can be written 
in terms of the scattering matrix at the Fermi energy. One arrives at the Landauer-Buttiker 
formula 



9ij = j^Tij, for i ^ j] gu = — [Tu - N] (2.4) 
where T^j is the total transmission coefficient from the ith to the jih lead, 

Tij = Tr \tijtlj^ ; 

tij is the scattering matrix between the states of the ith and the jth probe and is total 
number of scattering states (or total number of edge states in the presence case) at the Fermi- 
level. The above formula best illustrates the no-local aspect of quantum transport and has 
been instrumental to our understanding of Anderson localization and to the formulation of 
random matrix theory for quasi one-dimensional disorder systems. The DC conductance 
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Qij is well defined for finite-sized disordered region (embedded in an infinite open system), 
therefore are good candidates for scaling analysis; moreover, they are directly measurable 
in experiments. However, the formula is rarely used for microscopic calculations, because 
it introduces confinements and leads which can be cumbersome to address theoretically, 
particularly in the presence of magnetic field. Only recently it became understood that the 
presence of confinement potential at the edges changes the boundary condition for diffusion 
from {dn)(p = to {dn + ^dt)p = where t is the tangential and n is the normal direction 
of the edge [|13| . The boundary condition imposed by the perfect 2D leads will be discussed 
in this paper. 

In practice, the conductivity, a concept from classical physics, is still widely used, al- 
though for quantum mechanical system it no longer bears the local interpretation. It can 
be defined as the ratio between the average current density and the average field. Caution 
has to be exercised for finite-sized closed systems, where the DC dissipative conductivity 
is zero due to the discreteness of energy levels. The DC is usually asymptotically 
defined from AC conductivity Cxxi^) by taking the system size to infinite before taking the 
frequency (cu) of the external field to zero. Experimentally most of measurements of a^x 
and are deduced from longitudinal and Hall conductance(resistance) by extending the 
following classical relation to quantum mechanical systems: 

W 

9xx ^xx^^) 9xy ^xy-) (2*^) 

where g^x and gxy are combinations of the transmission coefficients. The "conductivities" 
thus defined are in fact conductances. The conductance or the transmission coefficients 
depend on a number of factors: the random scattering potential in the disordered region, 
the confinement potential, the property of the leads and even that of the thermal reservoirs. 
Result of a measurement is not independent of the way in which the measurement is done. 
Without checking that the measurement results are robust against size and geometry vari- 
ations, one can not be sure that the conductivity (or resistivity) are obtained. As we will 
show in section V, even in the perturbative regimes where quantum corrections are small, 
Rxx and Rxy exhibit the above Ohmic type of relations only under certain condition. For 
the critical regimes between the quantum Hall plateaus, experimental and numerical studies 
for different geometries have so far come up with a range of values for the critical value 
of axx HI], indicating the need for more careful study of finite-size scaling under different 



boundary conditions. We would like to address this point in our future study. 



III. LANDAUER-BUTTIKER FORMULA FOR RESISTANCE 

In this section, we express the longitudinal and Hall resistance of our system in terms of 
the total transmission coefficient. These Landauer-Buttiker type of expressions have been 
derived before and recently have been evoked by Shahar et al as a possible vehicle for 
the understanding of certain duality relation in the quantum Hall effect. Our derivation 
stress peculiarity of the measurement involved. 

The system as it is shown in figure 1 can not be used to extract both the Hall and 
longitudinal resistance. To do so while preserving the simplicity of the S-matrix, one can 
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attach additional voltage probes outside the scattering region. Suppose we attach 4 voltage 
probes and make a 6-probe Hall bar, as shown in figure 2a. The edge states outside the 
scattering region go from one probe to the next with probability 1 and all measurement 
result are related to one parameter, the transmission coefficient T = Tr{tt^}^ where t is the 
N X N transmission matrix of the right-going edge states over the disordered region. The 
transmission coefficient of the left-going edge states is also T by the unitary requirement of 
the S-matrix. The only non-zero elements among the T^jS are: T21 = T43 = T54 = Tie = ^5 
= Tes = T, Te2 = T^^ = where R = N -T. Passing a total current / in the 
x-direction, between probe 1 and 4, the resulting voltages can be solved (up to a constant 
shift) using the Landauer-Buttiker formula ( |2.4|) , bearing in mind that the total current in 
the voltage probes is fixed to be zero. The Hall resistance Rxy = ( V3 — V5)// = (V2 — Vq)/I in 
this particular setup is totally dictated by the property of the ideal leads since the voltage 
probes are across the leads. In the leads, the Fermi energy is pinned to its gap of bulk 
density of states. In such case, it is easy to show that the integrated current is equal to 
multiplied by the voltage difference of the two edges. Therefore, 

Rxy J^- (3.1) 

The longitudinal voltage drop can be shown to be 

]^^V^^ IR 

XX J INT 

We emphasize that while the above expressions hold for every realization of the disorder, 
they are not expected to hold for other geometries. As we have mentioned, the quantization 
of Rxy is a direct consequence of the measurement arrangement; it does not imply that Pxy 
is also quantized. Note that the quantization of Pxy implies that a^^ + a^y = Ncxy Were it 

to be true, this would give, for critical Gxy = cr^y = — 1/2, the critical axx = cr^^ = 
which contradicts with the current consensus that all integer quantum Hall transitions are 
the same and that a^^ should be independent of A^. This is one more reason that Rxx should 
not be taken as axx ioi N > 1. (This simple measurement setup is not able to detect whether 
N Landau levels are coupled by disorder or decoupled from one another in the scattering 
region, since it is only sensitive to the trace of the transmission matrix. Equilibrium among 
edges states are re-enforced by the thermal reservoirs attached to the probes.) 



IV. DUALITY RELATION FOR THE EDGE TRANSMISSION COEFFICIENT 

Recently, in studying of the critical transitions between the quantum Hall plateaus, 
Shahar et al 0,|1^ find that the longitudinal I-V curves near the critical magnetic field Be 
are non-linear and demonstrate a certain refiection symmetry with respect to the linear line 
at Be] more precisely, the longitudinal voltage and current appear to reverse roles at two 
filling fractions, u and z/^, on either sides of the quantum Hall transition: 

{VAu'),Uu')} = {Uu),VM}, (4.1) 
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with V and satisfying a definite relation suggestive of charge-flux symmetry in the bosonic 
view of the quantum Hall effect The authors point out that duality between charge 



and flux in the effective bosonic action can be the explanation of the observed relation ( [4.1|) , 
however there is one catch to this scenario- it requires the extra condition that the bosonic 
Hall resistivity remains zero across the phase transition, to which there has been no satisfying 
explanation. As an alternative the same authors also suggest a fermionic scenario, appealing 
to the Landauer-Buttiker expression similar to (|3.1| ) and ( |^ ) (their version is amendable 
to linear response). Noticing that for = 1, i?^;^; goes to its inverse as T ^ 1 — T = i?, 
Shahar et al propose that ( |4.1] ) can be explained within the Landauer-Buttiker framework 
if the following is true: 

T{E^ + A^;) = 1 - T{E^ - AE). (4.2) 

We show that the above relation can be the consequence of certain duality embedded in the 
Chalker-Coddington model for the quantum Hall effect. However, one can only write down 
the duality relation for the transmission coefficient of a symmetric sample with L = W. 



The Chalker-Coddington model [|TT[ consists of a lattice of directed links and scattering 
nodes (see flgure 3), representing the semi-classical orbits along the equipotential contours 
of the random potential and the tunneling among these orbits at the nodes. Each node is 
described by a 2 x 2 scattering matrix with random phases and a flxed probability Tq to 
scatter to the right and 1 — Tq to the left. Tq is a function of the Fermi-energy. In flgure 
3, we show one such network coupled to two edge states. At flxed Fermi energy, there 
are two equivalent ways to view the network: as one built out of clockwise guiding center 
orbits (the white squares) with probability 1 — Tq to tunnel to the neighboring orbits, or 
equivalently, of counter-clockwise orbits (black squares) with tunneling probability Tq. (At 
To = 1, the network breaks down to decoupled clockwise orbits and two edges state at the 
top and bottom; at Tq = 0, it breaks down to decoupled counter-clockwise orbits and two 
edge states at the left and right entrance of the sample.) The two states at energies related 
by T{E'f) = 1 — To{Ef) are dual in the sense that one ensemble can be mapped onto the 
other, if the system is inflnite, by reflection with respect to any of the discrete ridges in the 
direction of x±y. If the system is flnite and of arbitrary shape, the above symmetry is broken 
by the boundary. However, if the system is a square, two of the reflection axis survive. In 
this case the reflection along the diagonal brings the white squares onto the back square, 
and at the boundary, transmission channel becomes the reflection channel. Among many 
relations one can write down between the two states (or two phases), one is the following 

{T{l-To)) = {R{To)). (4.3) 

(In fact the distribution of T at 1 — Tq is identical to that of i? at Tq.) The self-dual point 
is apparently Tq = 1 — Tq = Tc = 1/2. Chalker and Coddington has made use of this fact in 
locating the critical point of the model. Linearizing the function Tq = T^^Ef) near Tc (Ec)^ 
E — Ec ^ (To — Tc)/Tq{Ec)^ the above gives equation (^?2|) and subsequently 



for (To — Tc)/Tc <C 1. (The above is rigorously true if the random potential distribution is 
symmetric with respect to positive and negative potentials. In that case it is appropriate 
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to assume Tq^Ec + AE) = 1 — To{Ec — AE). In this case there is no more need for the 
expansion.) If experiment is done in the geometry we consider, relation (|4.4|) wiU result as 
a simple consequence of the duality between the phase at Ef = Ef{To) and the phase at 

The experiment by Shahar et al was done in a Hall bar geometry with the voltage probes 
placed in the interior of the sample as shown in figure 2b. The S-matrix for the realistic Hall 
bar is significantly more complicated. Moreover the argument leading up to relation ( ^.4|) 
requires the refiection symmetry with respect to the two diagonal axis. There is no apparent 
reason why it should apply to a realistic Hall bar geometry with no such refiection symmetry. 
To find the possible connection between the experiment and the duality in the microscopic 
non-interacting models for the quantum Hall effect, one has to find out 1) whether or not 
a more general duality relation can be written down which is independent of geometry 
and boundary conditions 2) whether or not under certain conditions resistance measured in 
complicated geometries exhibit tractable scaling laws with some universal coefficients and 
exponents. 

The relation we wrote down for the discrete lattice model is somewhat artificial, since 
the network is an idealization of the mutually tunneling guiding orbits of arbitrary size and 
shape. It is then sensible to look for implication of duality in the continuous theories. The 
network model has been mapped onto the non-linear sigma model with a topological term 
|T5|] . The latter was shown by Pruisken et al to exhibit the appropriate asymptotic scaling 



property required of the quantum Hall transitions . If the network model is coarse-grained 
at length scale much larger than the lattice spacing, the clockwise and counter-clockwise 
orbits are lost, i. e. one can not tell the difference of the Tq and 1 — Tq state in the bulk, 
however differences does show up at the boundaries. Consequently, the bulk characteristic 
of the two phases, i, e. the diffusion constant or a^^ are the same, while their bare Hall 
conductivity a^y differ by one quanta. (It has been shown that the bare conductivity of 

the network model are a^^ = j^^^_^^~Jr^^,^2 and a^y = y2^(^°_yQ)2 P,P^). Since reffection is the 
operation that maps the Tq phase to the 1 — Tq phase, we consider how the non-linear sigma 
model transforms as under one reffection operation (x —x oi y ^ —y). The topological 
term, with a^y as its coefficient, changes sign while the a^^ term remains the same. Since the 
theory is periodic in a^^, the corresponding dual phases are that parameterized by a^y and 
n — a^y^ with n = . . . , —2, —1, 0, 1,2,.. .. (Note: a^y and n — a^y correspond to two filling 
factors or two Fermi-energies). Using the renormalization equations given by Pruisken et al 
0, one can verify the following relation between the renormalized axx^ (Jxy 



xy) ~ ^xx\'^ ^xy) 

^xy{(Txy) =n- Gxy{n - a° ), (4.5) 



^xx{(yly) =(^xx{n-a^^J 

'xyy^xy) = ^ ~ O-xyi"^ — O'xy 

if (Txxi^xy) = ^xxi^ ~ ^xy)- The above does not necessarily lead to the desired relation 

Pxx{(Txy) =Pxx{^-^Xv)' 



To obtain the above, one more constraint is required between axx and axy at the same a^y 
(filling factor, Fermi energy). For example, the above relation is satisfied if axy and axx obey 
(|4.5|), as well as the semi-circle relation 



7 



+ (^x. - 1/2)' = 1/4, (4.6) 

which is equivalent to 

Pxy = 1. (4.7) 

(Up to this pint, other form of the constraint is also possible. See later discussions on 
this issue.) Therefore, in both the fermionic and the bosonic picture for the quantum Hall 
effect duality alone is not sufficient to explain the experimental result ( ^TTI) . Both approaches 
require one more constraint: constant or vanishing Hall resistivity. As we explained before 
and will make more clear in the calculation to follow, the quantization of Hall resistance in 
our particular setup does not necessarily imply the quantization of Hall resistivity. It has 
been shown that equation (|46| ) or (|47|) hold for a classical version of the network model 
with no phase interference [|l6|,|l^,^. There are also numerical evidence that it holds for 



the original quantum version of the network model [|T8[. Theoretically the reason of the 
constraint is not clear. 



V. RESISTANCE IN TERMS OF THE BARE CONDUCTIVITY, 
PERTURBATIVE CONSIDERATIONS 

We next calculate the transmission coefficient, Rxx and Rxy in our geometry in terms of 
and axy in the perturbative limit. (The Chalker-Coddington model for one Landau level 
has very small bare conductivity a^^ < 1, therefore it can not be treated perturbatively.) 
We perform our calculation for high Landau levels (A^ > 1) and for the short-ranged random 
potential model, of which the bare conductivity a^^ can be large. The calculation serves to 
strengthen our view regarding duality in quantum Hall systems, it also has an interest of its 
own, since the problem of edge state transmission in the multi-scattering, diffusive regime 
has not been analyzed before. Our treatment is at most phenomenological. Its rationale is as 



following. In our previous work we have calculated the bilocal conductivity cr^zy(r, r^) to 
leading order in l/Cxx- that to leading order the non-local relation between current 

j(r) and external field E(r) can be decoupled and reduced to the familiar Drude equation 
(with some modification to be specified later) 

Uv)=a%£,iv), (5.1) 

where f = E + V/i/e is the electro motive field (/i is the local chemical potential). In the 
same work we also give the generating function from which one obtains the full quantum 
mechanical conductance. The effect of the edges can be treated within the non-linear sigma 
model, with a^^ and a^y as the coupling constants and the presence of the topological term 
alters the boundary condition for diffusion. In treating the quantum mechanical system as 
though it was classical, we are assuming that the field theory is renormalizable, i. e., all the 
quantum corrections to conductance can be accounted for by replacing the bare coupling 
constants with the renormalized reversion. The assumption so far has met no contradiction 
at perturbative level The new ingredient is the discussion on the boundary condition 
imposed by the 2D perfect leads. We show that the computation of transmission coefficient 
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and of current distribution reduce to solving the Laplace equation with identical boundary 
conditions. 



In reference [|T7[ it has been shown using the self-consistent Born approximation(SCBA) 
T9|] , which gives the leading order in l/cr^x^ ^^^^ bilocal conductivity tensor is of the 



following form: 



a^^r, r')) = la^S^, + + ag^yj S{r - r') 



^'A + ^iy%,'d,' + 4i%AS{y -w)- s{y)) 



^lA - ^i'y^uu'dl, - ail^HU^iy' -W)- 5{y'))\ d{r, r') 



xy 



+0{l/al,l/L), 



(5.2) 



where rf(r, r') is the diffusion propagator satisfying the — V^rf(r,r') = 5(r,r'), a^^, a^'y and 
a^y'^ are the SCB A version of the Streda conductivities p0|,|7|,p^ . The higher gradient terms 
are of higher order in //L, where / is the mean free path and L is the system length. The 
physics implied by the leading order expression is simple. The first term, the contact term, 
arises, in terms of the Drude's picture for conduction, from electrons accelerating in the 
combined external field E — ev x B. The second term is the diffusion term and it arises 
from the charge density fiuctuation Sn^ in response to the external potential, which is long- 
ranged. A local relation between the current and the eletromotive field can be obtained if one 
introduces an effective local chemical potential /x(r) to accounts for the density fiuctuation 
Tf\\ . However, equation does not quite recover equation 



T|) . One notices that in the 
diffusive term, a^y splits into two parts, a^^^, which appears in the bulk diffusion current and 
dxy^^ which is proportional to the edges current density and shows up only at the refiecting 
edges {axy = cr^y + ^i^)- It was shown that a^y is proportional to the rate at which the 
total density changes with the magnetic field at fixed Fermi energy |^ . In the low field limit 
^c^o <^ 1 (^c is the cyclotron radius and tq is the zero field elastic scattering rate), the 
Landau levels merge into a continuum and the density of states is only weakly dependent 
on the magnetic field, therefore the a^y term can be ignored. This is not the case in the 
high field limit with uJcTq > 1, when separate Landau bands are formed, and the density of 
states oscillates with the field. In general equation ( |0|) should be replaced by ^L7\ 



+ 



Wxx^lJ^^ 



^xy ^^^\ 



W) - 5{y)%^, 11/ e. 



(5.3) 



in which we have included an edge current /g = ajy' ji/e^ which is nothing but the extra 
edge current induced by an increase in chemical potential. 

In the disorder-free leads, = and = a^^y^ = is quantized, but is not 
always quantized in the disordered region 0]. It is apparent that when the edge current in 
the lead and the sample are different there has to be an edge current along the border with 
the leads amounting to = [cr^^y^ — a^y'^]iLi/e^ in order to satisfy current conservation at 
the corners. Therefore, for the geometry under consideration the local equation is further 
modified to be 



J, = alE^ + [alS,, + ai';e,.]d.^i/e + ai'/[S{y -W)- S{y)]S,,, fi/e 



lead 
xy 



(5.4) 
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To calculate the resistance, there are two approaches. One approach is to find the static 
electro-chemical potential and subsequently the current distribution, by requiring V x E = 
—d'Q/dt = 0, in which case, one can write E(r) = —V(f{r). The current conservation 
condition for the static case is V • j = 0, which gives in the bulk 

V2(-(^ + /i/e) = 0. (5.5) 



At the top and bottom edge, we have shown fT^] that current conservation leads to the 
following boundary condition 

aldyi-^ + fi/e) - al^d,{-^ + /i/e) = 0. (5.6) 

In other words, the eletromotive field E has to make an angle Q = tan~^ a^^/a^^ with the 
edge along x-direction: 

S (7^ 



XX 



At the left and right border with the leads, one can derive the following boundary condition 
by constructing a surface that encloses the edge current and imposing current conservation: 



e 

•lead 



3.-dyI'M=jT- (5-7) 

In the leads where a^x = 0, one can write down the following equation for the local current 
density [^]: 



Aead ^ ^lead^^^E^ ^ al^^^S^Mv - W) - S{y)]i,/e. (5.8) 



Combining (|5^,(l5^) dU), we get 



aldA-^ + ^/e) + - a':f^]dy{-^ + /i/e) = 0, (5.9) 

i. e., f has to make an angle 6' = tan~-^(cr^^ — cr^xy^)/cfxx with the border with the leads 
(along — y-direction ): 

^x ^ xy ^ xy 



^y ^xx 

Thus the problem of finding the resistance to leading order of l/cr^^, is that of solving the 
Laplacian equation for the electrical chemical potential with two tilted boundary conditions. 
Notice that as far as the conductance or resistance is concerned, one can arrive at the correct 
equations by using the simpler but incorrect local relation ( |5T|) . 

Alternatively, one can follow the Landauer-Buttiker approach, i. e., one finds the trans- 
mission coefficient. This also requires the computation of the bilocal-conductivity tensor. 
For our simple geometry, 

T = / dyrj dy2(jU^i,V2) (5.10) 

where Si and S2 are two the cross-sections at the left and right end of the scattering region. 
Making use of the SCBA expression ( |5.2|) for cr^iy(r, r^) , we get 
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1 nW rW 

T = --r dy dy\ald^ + aldy){alA-al&^)d{ry) 



\x=^,x'=L 



XX 

-^[d(L, W] 0, W) + ci(0, 0; L, 0) - d(L, VT; 0, 0) - d{L, 0; 0, W)]. (5.11) 

(J 



In reference [|T3, we have shown that the diffusion propagator satisfies boundary condition 
( |5.6|) at the reflecting edges, with d replacing —Lp + ji/e. One can easily check that it also 
satisfies boundary condition of the form (|5.9| ) at borders with the leads. The above expression 
can be simplified if one takes into account the additional detailed boundary condition that 
the returning probability via the incoming links is zero, the coarse-grained diffusive version 
of which is c/(L, 0; r) = rf(0, W] r) = 0, and the fact that transmission probability of left-going 
and right-going edge states are the same, i. e., (i(L, W] 0, Vl^) = rf(0, 0; L, 0). We get 



2^^^d{L,W-Q,W) (5.12) 



^XX 



Defining (/)(r) = 2^^2f^d(r; 0, H^), satisfies (pj), <^ and in addition the follow- 



XX 



ing initial condition (/)(L, 0) = 0, and 0(0, W) = since in the limit L ^ 0^ T ^ N. 
Thus within the diffusive limit the Landauer-Buttiker approach and the current distribution 
approach boil down to the same mathematical problem. Since the problem contains non- 
self-adjoint boundary conditions, its solution in the general case requires a special method, 
which we will discuss in the next section. However, some conclusions can be made before 
solving the equation. 

Again the fact that the Hall resistance is quantized can be demonstrated using 
boundary condition (|5^ ) alone. The current across the border with the leads is Ix = 

dy{a^^Sx + CTxy^y) the transverse drop in electro-chemical potential at the borders 
is A(j)y = dySy{x^y)x=o,L' Making use of the relation between Sx and Sy from ([0| ), the 
result Rxy = is recovered. 

The next conclusion follows from the two boundary conditions. Since the field S at 
the refiecting edges and at the borders with the leads is required to point at two specific 
directions, the field and current distribution is uniform only when the two directions coincide, 
i. e. when the bare conductivities satisfies the following constraint: 

(^L)' + <(^°,-^'^"') = 0, (5.13) 



which is equivalent to 



Pxy 



^ 1 



lead ' 
xy 



For a^^^^ = 1, the above gives the same semi-circle constraint ( |4.6|) of Ruzin et al. Under 



xy 



condition ( |5.13|) , the field distributions takes the simple form: 

cl)o{x,y) = -Eo{x + ^y), So{x,y) = Eo{x + ^y), (5.14) 



^XX 



where £"0 is a constant. The total current, the transverse and longitudinal electro-chemical 
potential differences can be easily calculated. Indeed Rxx.scba = ^Pxx- 
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The a^^ and a^y of the Chalker-Coddington model satisfies the above. However such 
bare values were computed by forcefully leaving out the quantum interference of the original 
model |T^. Since a^^ < 1, the model has no perturbative regime, i. e., there does not exist a 
length scale within which these bare values are the leading contributions. For the Gaussian 
white noise potential model a^^ = {2N — l)7r~^ sin^ a and a^y = N — 1 + a/n^ where a is 
a function of the Fermi energy (we do not need its expression here, see 0), they do not 
satisfy the constraint. In this case the field and current distribution is not uniform, Rxx is 
not proportional to Pxx and Rxy is not equal to Pxy 

As we mentioned before the geometry we consider (figure 2a) do not correspond to those 
used in the experiments, we would like to know how Rxx and Rxy change if we move the 
voltage probes into the scattering region as shown in figure 2b. Assuming that the voltage 
probes are very small in dimensions (much smaller than the sample length L and width W)^ 
we can treat the probes as perturbations in boundary conditions and get a rough picture by 
looking at the field distribution in the main strip in the absence of these probes. If condition 
( |5.13| ) is met and the field and current distribution is uniform, the electro-chemical potential 
drop transverse to the current is the same everywhere along the strip and the longitudinal 
potential difference is in proportion with the distance between the probes. Such systems 
have a quantized Hall resistivity, as it is the case with the Chalker-Coddington network with 



no phase coherence or a random version of it as discussed by . In the general case when 



the field is not uniform, the Hall resistance departs from the quantized value as we move 
the voltage probes to the interior of the sample. However if the Hall bar is long and narrow 
(L ^ W)^ the field distribution in the interior far from the ends takes the form of equation 
( |5.14| ), forced by boundary condition ( |5.6|) . If the probes are placed in the region where the 
field is nearly uniform, Rxx and Rxy are roughly proportional to the intrinsic p^^ and p^y. 
In this way the effect of the perfect leads is removed. 



VI. ANALYTICAL SOLUTION OF THE BOUNDARY PROBLEM 

To solve the boundary problem we use the method of conformal mapping, which was first 
used by Girvin and Rendell [^] to find the two-probe conductance in the case of absorbing 
leads. The idea is as following. Let us consider a parallelogram in the z^-plane {z^ = x' + iy')^ 
with its top and bottom edge parallel to x^-axis and its left and right side tilted from the 
vertical direction hj 56 = 6 — 9' -\- QQ. In this geometry, the uniform field distribution of 
( |5.14| ) satisfies both boundary conditions. The solution for the rectangle geometry in the 
z-plane (z = x + iy) can be found if one finds a conformal mapping z' ^ z which transforms 
the parallelogram to the rectangle. We start by writing down the complex potential (t){z') 
which can give the correct physical field ^o(^^2/0• The field are related to the complex 
potential through 

_ p , p 

dz ^ 

It is easy to see that the following complex potential 

(/)(/) = -lE^e-'^z'. (6.1) 

gives fo(x^y^) in the z'^-plane. The conformal transformation from z-plane to the z^-plane 
rotates the angle at x = 0, L by 59^ i., e., the rotation defined by 
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dz 

has the boundary conditions 

Im[f{z)] = 56, for x = Q, L; Im[f{z)] =0 for y = 0, W. 
The transformation has been worked out by Girvin and Rendell to be 

f{z)= J2 —{smh[nn{z-L/2)/W]/cosh{nnL/2W)}. 



n=odd 

The potential drops are: 

rL 



£"0 cos 



(9) / dxexp 
Jo 



— sinh[n7r(x - L/2)/W]/ co8h{n7rL/2W) 

.n=odd 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



/■W 

A(f)y = —Eosin{d — 56) / dyexp 

Jo 



456 



^ tanh.{mrL/2W)cos{mvy/W) 



nn 



The transmission coefficient is 



T = 4>{L, W) = N 



Ac/., + A<j)y 



(6.6) 



(6.7) 



In the hmit W = L, one can show that the two integrand in (|6.5|) and (|6.6|) are equal, the 
expression simphiies to be 



-L=W 



N 



sm 



56) 



N- 



cos{6') 



cos 6 + sin(6' - 56) cos{6') - cos 6 



(6. 



The resistance Rxx can be obtained by either plugging the above into the Landauer-Buttiker 
expression ( p. 2]) or by calculating the total current from the field distribution, then dividing 
A(px by it. 



VII. CONCLUSIONS AND DISCUSSIONS 



Prompted by the experimental evidence of duality near quantum Hall transitions, we 
probe the implication of duality in the Chalker-Coddington model within the fermionic non- 
interacting picture. For finite systems with boundaries, a duality relation holds for the edge 
state transmission coefficient, however, only for a class of systems with reffection symmetry. 
For such systems duality leads to measurable inverse relations of the longitudinal resistance of 
two states at opposite sides of a quantum Hall transition. In general situation (for arbitrary 
geometry) the duality translates to a diluted version cfxyi^xy) + ^xyi"^ — (^xy) = Oni 



key difference with previous authors on the subject [|T0| is that we think that the Landauer 
Buttiker argument constructed for resistance in symmetric geometry can not be expected 
to hold in other geometries, such as a Hall bar, since our calculations shows that resistance 
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is not necessarily resistivity. If experiment of Shahar et al is correct and the reversal of 
role of longitudinal voltage and current is indeed observed in the quantum scaling regime, it 
points to the possibility that the renormalized a^^ and a^y satisfy the constraint ( [4.6|) first 
suggested by the numerical work of Ruzin et al. This has so far not been understood within 
any theoretical framework. 

Our work unravels the mystery surrounding the Landauer-Buttiker argument. We per- 
formed a perturbative calculation for the edge state transmission coefficient in terms of the 
conductivities and analyze the current distribution. The presence of the 2D perfect leads 
gives rise to a tilted boundary condition for diffusion, which is similar and dual to the one 
at the reflecting edges, (dual in the sense that one boundary condition becomes the other 
as (j^y n — dxy)- Only when and satisfles the semi-circle constraint (equivalent 
to the quantization of Pxy)^ the current distribution is uniform and the resistance Rxx,xy 
is in proportion to the resistivity Pxx,xy Om phenomenological treatment is based on our 
knowledge of the bilocal conductivity cr^^(r,r') in the perturbative regime. To access the 
finite-size scaling of the resistance in the non-perturbative critical regime, one needs to study 
the form of cr^zy(r, r^) afresh. 

It is worth noting that we have considered the consequences of duality only within the 
zero-temperature, disordered and phase-coherent models for the quantum Hall transitions, 
while experiments are done at finite temperature where the electron-electron interaction (or 
the interactions between quasi-particles in the case of the fractional quantum Hall effect) 
must be considered. At a preliminary level interaction gives rise to a finite phase- coherent 
length which acts as a size cut-off in the renormalization of the conductivities when the 
elastic scatterings dominate over the phase-randomnizing inelastic scatterings. However, 
whether or not interactions are relevant at the quantum Hall transitions is an open question 
2M and there are indications that the non-interacting picture for the transitions do not em- 



brace all experimental facts (see ?? and the references therein). There have been interesting 
works which reveal non-trivial finite-temperature transport properties near critical transi- 
tions where interactions play key roles [^,0 . An explanation for the observed duality may 
lie in the understanding of the full problem with both interaction and disorder. 
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FIGURES 



FIG. 1. Edge states scattering through a disordered region. 

FIG. 2. The Hall bar geometry, a: fictitious version with the voltage probes set outside the 
disordered region, b: more realistic version with the voltage probes set in the interior of the sample. 

FIG. 3. Duality in the Chalker-Coddington network model for the quantum Hall effect. At 
fixed Fermi energy E'j, the system can be viewed as a network built of counter-clockwise orbits 
enclosing the shaded areas with tunneling probability 1 — To(£'/), or equivalently, a network built 
of clockwise orbits enclosing the white regions with tunneling probability To(Ef). The phase at 
To{Ef) can be mapped onto the phase at Tq(E'j^) = 1 — To(Ef) upon reflection with respect to one 
of the diagonal axis. 
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